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Abstract 

We investigate analytically the fermionic fluctuation determinant at finite temperatures 
in the minimal standard model, including all operators up to dimension 6 and all contri- 
butions to the effective potential to all orders in the high T expansion, to 1 loop. We 
apply the results to find corrections to the Sphaleron erasure rate in the broken phase. We 
conclude that the analytic treatment of fermions is very reliable, and that there is a great 
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deal of baryon erasure after the phase transition for any physical Higgs mass. 



^ ! 1 Introduction 

in 

It has been known since the pioneering work of Kirzhnits and Linde [Q and Dolan and 
Jackiw H that, at high temperature, the thermal populations of massive particles exert a 



symmetry restoring force on the Higgs condensate, so that at sufficiently high temperaure 
(T ~ mjj/g) the Higgs field loses its condensate and electroweak symmetry is restored. 
Under these circumstances baryon number is readily violated f||, |J. As temperature falls 



there is a symmetry breaking phase transition, which is first order for small ran- Presum- 



ably the universe underwent such a phase transition shortly after the Big Bang, and it is 



believed likely that the baryon asymmetry of the universe was generated at this time. 

Equilibrium thermodynamic information important to the physics of this epoch, such 
as the effective potential and the free energy of certain saddlepoint solutions, can be com- 
puted in the Matsubara (imaginary time) formulism. The particular calculation we will 
be interested in is the Sphaleron rate in the broken phase shortly after the phase transi- 
tion. This was first investigated at tree level in [E], and one loop corrections have been 
computed within a high temperature approximation in || |7] . This approximation assumes 
that the contributions of nonzero Matsubara frequencies can be absorbed into shifts in the 
parameters of the zero Matsubara frequency modes, which then constitute a "dimension- 
ally reduced" three dimensional theory ||. The same idea forms the basis for an extensive 
investigation of the effective potential and the strength of the phase transition being carried 
out by Farakos, Kajantie, Laine, Rummukainen, and Shaposhnikov (FKLRS) || [10 . 



Why is it necessary to include 1 loop corrections? Normally, one loop corrections modify 
tree level results by only a few percent, but there are exceptions. The finite temperature, 1 
loop correction to the Higgs mass is significant because, although the correction is paramet- 
rically order g 2 , the mass is a dimensionful parameter, and the importance of the thermal 
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effect is enhanced by T 2 /m 2 , which can be large. The phase transition temperature is the 
point where the thermal, 1 loop effect and the tree level effect are almost equal. Finite tem- 
perature calculations also have potential infrared divergences from loops containing zero 
Matsubara frequency bosons, and it is not clear that perturbation theory always works; 
the object of the dimensional reduction program is to separate this problem from those 
effects which can be treated perturbatively more reliably. Finally, loop corrections can also 
be important when a coupling is unnaturally small; for instance, the Higgs self coupling A 
recieves a 1-loop correction which goes as gf (g t is the Yukawa coupling of the top quark. 
In this paper we will always take g t = 1, corresponding to m t ~ 174GeV). This correction 
does not depend on A, and if A happens to be small then the correction can be very im- 
portant. This is the case in the standard model whith a light Higgs boson and a heavy top 
quark. It is quite necessary, then, to consistently use 1 loop relationships when describing 
anything directly related to A, such as the physical Higgs mass and the (finite and zero 
temperature) effective potential. In this case, however, two loop effects are still expected 
to make a small perturbative correction to the one loop effects, so the one loop calculation 
should be reliable. 

With this in mind, Diakonov, Polyakov, Sieber, Schaldach, and Goeke (DPSSG) have 
made a direct numerical evaluation of the fermionic fluctuation determinant and have 
concluded on its basis that, when the top quark is heavy, fermion fluctuations make the 
Sphaleron rate very different from the high temperature estimate |TTJ. This brings the 
dimensional reduction program into question, and bears further investigation. Here we 
give a one loop, perturbative treatment of the fermions in the presence of zero Matsubara 
frequency background fields. The free energy of the fermions can naturally be expressed in 
terms of operators made up of zero Matsubara frequency bosonic fields and their deriva- 
tives. In section II we compute these, including induced masses (dimension 2), corrections 
to couplings and wave function renormalizations (dimension 4), and nonrenormalizeable 
dimension 6 operators. The expansion in operator dimension appears to be very well be- 
haved when the background bosonic fields are slowly varying. In section III we apply these 
results to a computation of the Sphaleron energy. In section IV we use this calculation, 
along with 1-loop relationships between the couplings and the physical Higgs mass and a 
calculation of the bubble nucleation rate which determines the temperature at which the 
phase transition occurs, to investigate the erasure of baryons after the phase transition is 
complete. In the last section we draw conclusions. Some technical details are relegated to 
an appendix. 

2 Integration over fermions 

We will work in SU(2) Higgs theory coupled to all the fermions of the standard model. 
Neglecting hypercharge significantly simplifies the picture whithout profoundly changing 
it, because the Weinberg angle is small and because the fields which make up the Sphaleron 
are almost exclusively weak isospin, and not hypercharge, fields. (It is known that, at tree 
level, including hypercharge with the physical Weinberg angle only changes the Sphaleron 
energy by about 1% ||). The Lagrangian, suppressing summations on generations and 
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colors, is 



C = £ h + £ f 

C b = ±F^ + (D^D^-ml&<S> + \(&<S>) 2 

(We use a Euclidean metric g^ u = 5^ u and Euclidean 7 matricies which satisfy the algebra 
7 M 7„ + 7i,7 M = 25fj, u .) The fermion kinetic terms apply for every doublet and singlet, but 
the mass term only applies for the top quark. Actually the mass term as written is for a 
bottom, not top, type quark. This is for notational simplicity only; to make the top quark 
massive one should systematically replace $ with ir 2 &* in all that follows; our conclusions 
are completely unchanged. 

In the 1 loop approximation we ignore the fermions' coupling to nonzero Matsubara 
frequency excitations of the bosonic fields. The path integral over the fermions is gaussian; 
its contribution to the partition function is DetH, and the corresponding contribution to 
the effective action is — Tr \nH, where H, written as a matrix acting on [■?/>£ i/)r] t , is 



H 



9t& a 



7^<W 
7/^ 



+ 



9t& a 



H + H T , (1) 



where a and (3 are SU(2) indicies. 

Our analysis will be based on an expansion of — Trlnif in Hj. To illustrate the idea 
of expanding the log, consider a simplified example in which the gauge field is everywhere 
zero. In this case we have a Dirac fermion with a (spatially varying) mass, H = 7 M <9 M 



m, 



where m 2 = g 2 &<&. We assume that m 2 is smaller than the lowest eigenvalue of — d 2 , in 
which case it is legitimate to expand the log. This will generally be the the lowest 

eigenvalue of — d 2 is set by the square of the lowest possible Matsubara frequency, which 
for a fermion is (nT) 2 . The log becomes 



Tr ln(7^ - m) = -Tr In 7^ - ]T ^ Tr 



n=l 



n 



m 



(2) 



where as usual 1/(^8^) is defined as 7 At <9^/<9 2 . We get Feynman diagrams by the usual 
trick of inserting complete sets of states and Fourier transforming. The first term in the 
sum is a divergent vacuum energy and should be removed. All terms with odd powers of m 
vanish when we take the trace on Dirac indicies. When the mass is position independent, 
the resulting terms are 



rT- 1 p 00 

J dx J d 3 x ^2 (■ 



4m 



2n 



2n 



ko 



d 3 k 



(2tt) 3 (k 2 + k 2 y 



(3) 



We have used the shorthand Y,k ^° mean a sum in which k takes on odd integer multiples 

of 7rT; we will use this shorthand throughout. Also, k 2 means k 2 and k means The 
(— ) n arises because when we Fourier transform d — > ik. The 4 is from the Dirac trace, and 
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2n is the symmetry factor of the diagram; we saw how it arises in the expansion of the log 
above. To get the free energy density we drop the space integral. 

The first two integrals in the series are ultraviolet divergent and must be performed with 
some care. We have done so in the MS scheme by first performing the sum on Matsubara 
modes and conducting the spatial integrals in 3 — 2e dimensions. The higher order terms are 
convergent and may be performed directly by doing the integral over d 3 k first. Performing 
the integrals, the free energy per unit volume is 

2 T2 4 2 7i; -21n 7 r + ln( ; u 2 /T 2 ) » 4m 2 " (1 - 2 3 - 2 ")C(2n - 3)I> - f) 
m 12 16vr 2 ^ ' 2n 4vr 2 «- 2 T 2 ™- 4 r(i)r(n) 1 ] 

with the usual MS renormalization point. This is the 1 loop fermionic contribution to the 
finite temperature effective potential, to all orders in the high temperature expansion^. Note 
that, as expected, the sum is a Taylor series in m 2 with radius of convergence (vrT) 2 . By 
m t we mean the thermal top quark mass given by m 2 (T) = g 2 $^$(T), not the vacuum top 
quark mass m 2 (0) = g 2 &Q(T = 0). Near the phase transition temperature the difference 
is quite substantial. In the simplest approximation $^$(T) = $t$(T = 0)(1 — T 2 /T 2 ), but 
near T c it is necessary to treat the influence of infrared bosons more carefully. As we will 



see in Sec. IV, at the temperature of interest the Higgs VEV u = y$t$/2 < T and for 
g t = 1 we find (m/irT) 2 < l/(27r 2 ), so the convergence of the series is excellent. Of course, 
well below the phase transition temperature it is no longer true that m t « nT and our 
approximation scheme will fail. However, at these temperatures the Sphaleron energy is so 
large that the baryon erasure rate is utterly negligible, so our approximation scheme should 
apply in the interesting range of temperatures. 

Another way of understanding the convergence of the power series is to think of the 
different Matsubara frequency contributions of the top quarks as distinct, very massive 
species of a three dimensional theory, with masses given by the Matsubara frequencies 
m 2 M = ((2n + l)7rT) 2 ; we are then expanding in the ratio of g 2 &$> (and eventually in 
products of other infrared bosonic fields like A { , and their derivatives) to m 2 M . Such an 



expansion has been considered in a theory with only fermions and a scalar in |12|, |13 
where it is also concluded that such an expansion is very accurate. An expansion like the 
one used here is not justified for bosons because symmetric boundary conditions in time 
make the lowest eigenvalue of the operator d 2 is zero; or equivalently the lowest Matsubara 
frequency is zero, so we cannot expand in the ratio of a field value to this Matsubara 
frequency. (However we could expand all the other Matsubara frequencies in the way 
described here, see [f|.) 

For simplicity of notation, in the remainder of the paper we write 

T ^ f d 3 k 1 _ 2 7g -21n7r + ln(/x 2 /r 2 ) _ 

^ I (2vr) 3 (k 2 + kl) 2 16tt 2 " 1 ' 



fe 



and 



T yflL 1 ^ (l-2^ K (2iV-3)r(iV-j) 



2 An integral expression for the fermionic contribution was first found by Dolan and Jackiw g , who also 
found the first two terms in the expansion given here. What we have done is found the complete Taylor 
series for the known integral expression. 
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in particular, 

When the mass is spatially varying the first term in the sum in Eq.(^) becomes 
f -1 f d 3 p f d 3 k tr(7"7")fc„(p + fc)„ 

If we assume that p 2 < k 2 + k" 2 , which means that the mass varies on a scale large compared 
to 1/ 7rT, then we may expand the denominator in a geometric series and extract a power 
series in p 2 . We can then Fourier transform to position space to express the result as a 
derivative expansion for m. The resulting free energy is 

£ x (!l m i + DA (Vm) 2 - — (V 2 m) 2 + -(VV 2 m) 2 ..) . (9) 

y 12 3 10 J 

In realistic cases the Fourier transform of m lies primarily at small values of p but has a 
rapidly decaying exponential tail which goes above p = ttT. In this case the series will 
be asymptotic and its reliability will depend on how rapidly the exponential tail falls off. 
For the Sphaleron at T ~ T c we expect the exponential tail of the Fourier transform of 
the gauge and Higgs field configurations to fall roughly as exp(— p/mw(T)); at T ~ T c , 
mw{T) = g w v{T)/2 is << nT and the convergence of the (asymptotic) derivative expansion 
should be excellent. Of course we will check this by explicit calculation. 
For the more general H the lowest nonvanishing term will be 

1 #z Q Hj_ 1 rjh_ T y f_iL ( ^ m ik,i{k+p)MYl u ) 

2 iT Hf H 2Ho ^2J (27r)3 J ^i (27r)3^ W P)) (k 2 + k 2 )((k+pf + k 2 ) 

+ 2i (2703^/(2^ [ 2 ^(P)^ «J (F + ^)(( fe+p ) 2 + ^) ^ 

which corresponds to the Feynman diagrams illustrated in Fig. III. Higher order contribu- 
tions can be gotten similarly by going to higher powers in Hj. 

The contribution of fermionic fluctuations to dimension two and four operators in 
SU{2) x U(l) Higgs theory in MS with realistic couplings has recently been worked out 
by FKLRS |T0|| . We independently performed the calculation in SU(2) Higgs theory; our 
results concur. We have also extended the calculation to find all nonrenormalizeable oper- 
ators induced by fermions at dimension 6. Our results follow. 

The dimension two operators induced by fermions change the masses of the Higgs and 
Aq fields. The induced terms are 

N c g 2 T 2 t N d g 2 ,T 2 a9 

-S!± — $t $ + ^^ — A 2 (11) 
12 24 v ; 

The first term is largely responsible for the restoration of symmetry at high temperature. 
The second is the familiar Debeye screening mass. (Here and throughout iV c = 3 is the 
number of colors and Nd = 12 is the number of left handed fermion doublets.) 
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At dimension four, fermions introduce the following corrections to couplings and wave 
functions: 

N n 2 a 2 1 

N c g$DA ($t$) 2 + N c g 2 t DA (A$) t A$ + -^p(D4 - ^)K&® 

+^^F«F« + ^f(D4 - ^ 2 mA n Di A y - ^{Alf ■ (12) 

The coefficients of D4 give the fermionic contributions to the one loop beta functions and 
anomalous dimensions for the bosonic fields in SU(2) Higgs theory, and agree with well 
known results. If the dimensional reduction scheme is justified then these terms should give 
an accurate approximation of the fluctuation determinant. To check this it is necessary to 
go one order higher and see if the contributions of dimension 6 operators are small. 
At dimension six, fermions contribute the following terms to the effective potential, 

_ 2N^m^ f + N c gfgjD6 {Al)2 ^ ^ 

the following Higgs derivative terms, 

- ^($t$)(A$)t A $ _ ^|^ ($t A $)( A$ )t $ 

_^D6^ ($t$) ^ ($t$) _ ^D6 (jD 2 $)tjD2<|) ; (14) 
3 3 

the following mixed derivative terms, 

Nc9l f D6 Al(D^yD^ - N <&& m ' (Ai4o) , (Ai4o) ag t g , (is) 
the following A derivative terms, 

^^mlMAt) - —^^—(D 2 A ) a (D 2 A ) a , (16) 
the following mixed terms, 

8 21 
and the following gauge field terms, 



N c9t g w D6 ^ N d9w D6 2 N d g w D6 b b 



- - ^f|^(A^) a (AF fe ,) a . (18) 

Fermionic contributions to other independent gauge invariant dimension 6 operators 
(such as (Aq) 3 ) vanish. 

Exactly the same terms arise whether the massive quark is top type or bottom type. If 
both types are massive (in one generation) then each occurrence of g™ becomes g™ + g™ in 
the above, and exactly one mixed term appears at dimension 6, 

- AN ^ t9 * DQ &D j iT 2 $*{D j iT 2 $*y$ . (19) 
3 



Fermions also induce pure QCD operators and operators containing both weak and 
strong fields; these are unimportant here and are listed in Appendix B. 

We should comment that, except when they co incidentally vanish, the fermionic con- 
tributions to dimension six operators tend to be much larger than the contributions from 
the nonzero Matsubara frequencies of boson fields, some of which have been worked out 
in [ 14| . This is partly because the top quark is heavier than any of the bosonic degrees 
of freedom and partly because its lowest nonzero Matsubara frequency is 7rT, rather than 
2nT . Because of this the bosonic equivalent of the coefficient D6 is smaller by a factor of 
7. Hence, we anticipate that if the expansion in high dimension operators is well behaved 
for fermions then it will also be well behaved for bosons. 

We have also computed the fermionic contributions to masses, couplings, and wave 
function renormalizations in a slight modification of the proper time technique of DPSSG. 
The procedure is outlined in Appendix A. The results turn out to be identical to those in 
MS except that DA is modified to 

DA - 7*-21nvr + huy/r 2 ) , , 

^ ^proper time -^g 2 ' \ I 

where the definition of the scale /i is explained in Appendix A. This expression relates the 
MS and proper time renormalization points. The appendix also presents the calculation of 
the vacuum effective potential and the physical Higgs mass in the proper time scheme. 



3 The Sphaleron 

We want to apply these results to find the free energy of a nontrivial field configuration 
which solves the classical equations of motion. Klinkhammer and Manton have shown || 
that the classical equations of motion of the SU(2) Higgs system can be solved by an Ansatz 
of form 

x + iy 



v h(r) 

$ = — 

V2 r 



z 



A{ = ^-r k e ijk . (21) 



The lower index on A is the Lorentz index and the upper index is the group index, and h and 
/ are functions of r alone which are to be determined by minimizing the configuration's 
energy, subject to the boundary conditions /(0) = h(0) = and /(oo) = h(oo) = 1. 
(Here and throughout v is the Higgs VEV, v 2 {T) = 2$ t $(T).) This solution is called 
the Sphaleron, and when Es p h » T, the formalism of Langer tells us that most baryon 
number violating events should occur because of phase space paths which pass close to this 
configuration. Arnold and McLerran have applied this idea to estimate the rate of baryon 
number violation in the broken electroweak phase to be M 

(1Xu -1SN F T ^- ( %\ N tv NV TOt KeM-E Sph /T) (22) 



N B dt \47r/ 2m w \g w T 

with Np = 3 the number of generations. The prefactors arise from the zero and unstable 
modes of the Sphaleron, and are evaluated in [[15], |fl. At 1 loop k is a product of fluctuation 
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determinants around the configuration, 



_ Detff f DetKo Y 
K ~ Det/fo VDetK 7 / ' { 6) 

H is Eq. ([I]) in the Sphaleron background and H is Eq. ([I]) in the naive vacuum. K is 
the bosonic fluctuation determinant in vacuum, and K' is the bosonic determinant in the 
Sphaleron background, but with the zero and unstable modes removed. 

We should comment that the division between In k and —E sp h/T is somewhat arbitrary. 
In particular it is renormalization point dependent. We see this explicitly from our calcu- 
lation of the contributions from fermions to dimension 4 operators (Eq. (0)), which are 
part of — In k, but which have with coefficients oc D4 which explicitly depend on //. The 
calculation of Es p h correspondingly depends on coupling constants g w , A which depend on 
//. 

To compute Eg p h — TIuk we need to evaluate the tree level Lagrangian terms 
(<frt<3>) 2 , and F°jF°j, and the operators found in Section II, in the Sphaleron background, 
Eq. (pip. First note that Aq = 0, so all terms including Aq vanish. Effective potential 
terms give 

v 2 

$t$ = _/ l 2 ; (24) 
($t$) 2 = (25) 

($t$)3 = i 6 ... . (26) 
8 

We will also need the term ($t$) 3 / 2 = v 3 h 3 /(2\/2) which arises from bosonic fluctuations, 
and we should add a independent constant to make the global minimum of V(v) zero, 
to subtract out the energy density in the absence of a Sphaleron. 
The dimension 4 derivative terms are || 

f'2 f2f X _ f\2 

papa = 1 6 Z_ + 3 2 J \ A J) , (27) 



(A$) f A$ = y (2 k2{1 r2 /)2 + h> 2 
We calculate that the dimension 6 derivative terms are 



(2i 



UcF-F^ = 96 f2fi ^ f) , (29) 

f"2 f"f(l - f) f 2 (l - f) 2 

(D t F l3 r(D k F kj y = 8^-32 J Jl r4 J > + 32 J { ^ 1 > , (30) 

2 h 2 ( fl2 f2(i _ f\2\ 

&*F%F% = -^^ + 32 ;i r4 ^ j, (31) 

2 / A' A/l -f\2\ 2 



(DH>yDH> = V -(h" + 2^-2 fe(1 ^ ) , (32) 
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a i ($ t $)^($ t $) 



v 4 ( h 2 (l-f) 2 , /2 N 



v 4 h 2 U 2 
v 4 (h 4 {\ 

1 



(x 2 + y 2 ) + h 2 h' 2 



(33) 
(34) 
(35) 



Because / oc r 2 and h oc r at small r, all of these terms are nonsingular at r = 0. 
The last term, and only the last term, is not spherically symmetric. The departure from 
spherical symmetry arises because the Dirac equation in the presence of the Sphaleron is 



not spherically symmetric when only one flavor of quark has a mass [II]. DPSSG evaded 
this problem by giving both flavors equal masses, which restores the symmetry. As we have 
seen, giving both quark flavors masses introduces a new dimension 6 operator, whose free 
energy density in the Sphaleron background is 



u 4 h\l-ff 



(r 2 + z 2 



(36) 

g 2 , the z 2 combines with 



The coefficients of this term and ( |35| ) are such that, when g\ 
the x 2 + y 2 to restore the spherical symmetry of the terms. 

The approximation of DPSSG differs from the results with only the top quark massive 
by the contribution of Eq.(p6|) and is accurate only when this term is small. This is the 
case only when derivative dimension 6 operators give very small contributions, which is 
precisely the case where dimensional reduction is accurate. 

To compute the Sphaleron energy it is convenient to follow || and introduce a dimen- 
sionless radial distance £ = g w vr. The contribution to the Sphaleron energy from the 
effective potential is then 

^r«-^. (37) 

9w Jo g 2 w v 4 

The contribution from kinetic energy terms, including the fermions' contribution, is 




+ (l + JV^£>4) [h 2 {l-ff 



where derivatives are with respect to £. 
The contribution from Eq. (18) is 



(3? 



AnujO^gty 
g w 128vr 4 T 2 



180 



15 



96 



n2 



fV(l - /) 



^-32 



r/(i - f) 



32 



f\l-ff 



(39) 



We can get the other dimension 6 operators from Eqs. (|I4] - [18]) and Eqs. (^9| - |35|) by 
always replacing D6 with 7C(3)g^z/ 2 /1287r 4 T 2 , removing g w 's, replacing g t with gt/g w and 
r with £, and integrating over (47rv/g w ) J £ 2 c?£. For v ~ T, D6 is very small and dimension 
6 operators have the parametric appearance of two loop effects. 
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From the discussion after the introduction in Q we see that we want to find the config- 
uration of form ( |2TD with minimum free energy, that is the configuration which maximizes 
(Ilzero mode cont.)«exp(— Es p h/T). This is a formidable task, since we cannot compute all 
of these terms analytically; in particular the zero mode contributions and part of the zero 
Matsubara frequency bosonic contribution to k are beyond our analytic abilities. Fortu- 
nately, the Sphaleron is a saddlepoint configuration, and we should get almost exactly the 
right free energy if we include the dominant effects in the computation of the field configu- 
ration /(£) and h(£), and then compute the other corrections holding / and h fixed. This 
is because, as a saddlepoint, the free energy of the Sphaleron only changes quadratically 
with small changes to h and /. 

To illustrate this point consider the Sphaleron configuration computed from a tree 
level Lagrangian F^F^/i+X^Q — u 2 /2) 2 and suppose that we are only interested in the 
correction fermions induce in the gauge fields, A^g 2 DAF^F^/Yl. We will choose \i = 
mw{T = 0) in the proper time scheme and T = 100 GeV, so D4 ~ —0.0135. Solving for 
the Sphaleron using the tree action we find Es p h — 35.0975z/, and estimating the fermionic 
correction as 

£ f d( (> + S f \J )2 ) c -0.4433, (40) 

we get a total free energy of 34.6542z/. When we solve for the Sphaleron configuration, 
including the fermionic contribution as well as the tree terms, we find the Sphaleron energy, 
which now includes the fermionic correction, is 34.6520za The correction from this fermionic 
contribution is about 1% of Es p h, and the error in estimating it at fixed / and h is about 
0.01%, which is quadratic, as expected. 

Obviously, though, this will not do when a correction is actually substantial and the 
modification of the free energy is comparable to Eg p h- This is potentially the case for the 
fermionic correction to the (<E>t$) 2 term in the effective potential. For the above parameters, 
the coefficient of the correction is N c gfD4 = —0.040, which is to be compared with A(/x) = 
0.050. Because of an unfortunate choice for /a, the "correction" for top quarks is almost as 
large as the tree level term itself, and we certainly cannot trust a result in which / and h are 
computed with the tree term and used to compute the quark fluctuation determinant. The 
only consistent, renormalization point independent thing to do is to include those fermionic 
contributions which correct operators appearing in the tree level action in the calculation 
of the Sphaleron configuration, that is to use the ii independent quantity A(/i) +N c gfD4 in 
the calculation of /, h, and Es p h- This is particularly important for the scalar self-coupling, 
because as mentioned earlier it is unprotected from large radiative corrections and in fact 
the top quark contribution here is substantial. 

To work in terms of physical quantities, we first find a relation between the couplings 
at the scale \i and vacuum, physical masses. That is, we find \{p) in terms of m#(T = 0). 
We perform this calculation, including all one loop, fermionic contributions, in Appendix 
A. We ignore the bosonic corrections because we are most interested in understanding 
the fermionic radiative corrections in this paper, and because the bosonic corrections are 
smaller by a factor of 9m^y/12mf ~ 1/30. From the Appendix A results we find A(/i) = 
{m 2 H /2i>l) — (/V c ^/16vT 2 )(ln(/i 2 /m 2 ) — a), where a = in MS and a = j E in proper time 
regulation. The quantity we should use in calculating the Sphaleron configuration at a 
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temperature T is 

A T = A(//)+iV c <7 t 4 £>4, (41) 

which has no fi or renormalization scheme dependence, as we can note by inspecting Eq. 
(|) or Eq. ©0. 

A particularly convenient choice for /i is 

\i = exp(— je)^T, MS or \i = exp(— r y E /2)7iT, proper time (42) 

because in this case At = A(/x) and the Higgs field and gauge field wave functions take their 
tree values; but there is no need to choose this scale, as long as we compute the Sphaleron 
configuration and energy using couplings and wave functions corrected by the fermion 



contributions as in equations (|3qj4l[ ). This is precisely the prescription of the dimensional 
reduction program. We can test its reliability by seeing how large the remaining corrections, 
those coming from dimension 6 operators, are. 

Of course, it is also necessary to include corrections from zero Matsubara frequency 
bosons, which are expected to be considerable and cannot be evaluated with a derivative 
expansion, as we have already discussed. The largest part of this correction comes from the 
cubic effective potential terms. The residual correction when this is completely removed 
is In K ~ 1.5 0, |16|] which is small enough that we can trust the computation in terms 
of fixed / and h as discussed above; but the effective potential contribution is large (at 
T ~ T c it changes the very nature of the phase transition and should be considered an 
order 1 correction to Egph) and should be included in the evaluation of the Sphaleron 



configuration, as advocated in fL7 |. But before we can compute the Sphaleron energy 
including these corrections we must discuss the nature of the phase transition to find out 
at what temperature the baryon erasure begins. 



4 Phase transition, Sphaleron energy 

Let us briefly review the electroweak phase transition. At one loop the zero Matsubara 
frequency bosonic excitations generate negative cubic terms in the effective potential, which 
becomes (Qw — 0) |0 



v(v) = _9^H , (_< , (4^ 2 + 3^ + 8A)T 2 \ y2 _ ?T _ ( IIT 2 v*\ ^ 
1 ' 90 V 2 32 J 16vr 4tt \ 6 4 ) 

— h at— + dimension six , (43) 

127T 4 

where m\ = \rpu 2 /2 + V"(u = 0), m\ = Aj-z/ 2 + m|, g* = 106.75 is the number of radiative 
degrees of freedom, and "dimension six" means the higher terms found in Eq. (TJ) . At high 

3 The alert reader may notice that A(/u) also contains a term (m? H /2v%)(N c gt /167r 2 )(ln(yU 2 /m( ) — a), 
which apparently spoils the cancellation of the [i dependence discussed here. This is true, and it gets 
another correction because v Q is /i dependent; but the mass squared parameter TOq of the effective potential 
contains a similar /i dependence, so the effect of the correction is just to shift the location of the minimum 
of V by a proportional amount, without changing its height. The Higgs field wave function has the same 
dependence, so this has no influence on the Sphaleron energy. 
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temperature V has only one minimum at u = 0, but as temperature drops the negative cubic 
terms generate a second "asymmetric" minimum. At some temperature the second min- 
imum becomes more thermodynamically favorable, and bubbles of the asymmetric phase 
begin to nucleate and grow shortly thereafter. The temperature at which the nucleations 
become common can be computed by standard techniques p^| . The expanding bubbles 
liberate latent heat, so that the temperature after the transition is somewhat higher than 
the nucleation temperature. The temperature the plasma reheats to due to this latent heat 
is determined by the condition that the broken phase energy density E = V — TdV/ dT 
must equal the symmetric phase energy density at the temperature where the nucleations 
occurred. It is at this reheat temperature, immediately after the phase transition, that 
quasi-equilibrium erasure of any baryon number excess generated at the phase transition 
begins. It continues for all times thereafter, but as Eq. ( p2|) shows, the rate depends 
strongly on the Sphaleron energy, which changes rapidly with T, as v moves towards its 
zero temperature value. Almost all the baryon erasure takes place within a fraction of 
a Hubble time, so the Sphaleron rate is only relevant at temperatures quite close to the 
reheat temperature. 

We have computed the reheat temperature and the Sphaleron energy at the reheat 
temperature for a number of physical Higgs masses, using the 1 loop relations between 
Higgs mass and A presented in Appendix A and including the negative cubic effective po- 
tential terms from zero Matsubara frequency bosonic modes. (We always use g w = 0.65 
and take the Weinberg angle to be zero. In fact, to account for the quark contribution to 
the gauge field wave function, the value of g w we should use is g w ([i) — A^g^ZM/G, which is 
T dependent. The T dependence is very mild, behaving as g w * (A^g^,/487r 2 ~ 0.01) InT. 
Since this dependence is so weak, and since we are not including the influence of bosons' 
nonzero Matsubara modes, which will contribute an opposite and slightly larger tempera- 
ture dependence f9j, we will not worry about it. However we will take full consideration of 
the correction to A, which as discussed is not at all weak.) The results, together with the 
contributions from dimension 6 operators, are presented in table 1. The z/ 6 contribution to 
the effective potential was computed by inclusion in the calculation of the Sphaleron con- 
figuration, and the others were performed perturbatively. It is clear from the table that, 
as expected, dimension 6 operators make only a tiny contribution to the Sphaleron energy. 
Hence we can conclude that near the phase transition temperature the expansion is very 
well behaved and the use of the dimensionally reduced theory is well justified. The largest 
dimension 6 correction comes from the effective potential term, owing to the high power 
of g t 'i this and only this term may not be completely negligible, increasing the Sphaleron 
energy and the strength of the phase transition by a few percent for light Higgs. Generally 
the terms become progressively less important as they contain more derivatives. In partic- 
ular the very small contribution from (D 2 Q)^)D 2 Q gives confidence that the expansion of 
Eq. (H) in derivatives is justified. 

We have not continued the table down below 30 GeV partly because this range is 
experimentally excluded, partly because questions of vacuum stability become increasingly 
hard to avoid in this range, and partly because the results barely differ from those at 
itih = 30 GeV. This is because the one loop relations for the Higgs mass give nonzero At 
even as the physical Higgs mass ran — ► 0. 

The next step is to use these results to determine the baryon number depletion. From 
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Eq.(p2|) we find that baryon number is depleted by a factor of 

exp ( [ 13N F T 4 ^—N tr NV rot ( ^ V k exp ^^dt) (44) 
yt ieheat V4tt/ 2m w \g w T J T J 

The elapsed time is related to the change in temperature by dt = dT/(HT), with H 
the Hubble constant, which is approximately H = T( y 8% 3 g^/90) 1 ^ 2 (T/m p i). To perform 
the integral we must numerically repeat the evaluation of the Sphaleron energy at many 
values of T close to the reheat temperature. Most of the coefficients in Eq. fl4"3|) are 
given in [15|]. We use the values found there for uj_ and the M, even though they were 



computed for slightly different / and h, because the product of these factors turns out to 
be very insensitive to changes in the configuration ||15|| . The other factor we require is the 



difference between the Sphaleron energy which we have calculated and the (1 loop) value 
of Esph — Tin ac. This is due to derivative corrections from zero Matsubara frequencies and 
can be approximated from the results of Baacke and Junker [7|], who find that, when the 
full tadpole is removed (ie. all effective potential contributions are subtracted out), hiK is 
about 1.5, independent of v and quite weakly dependent on X T . (Again, this quantity was 
computed for different /, h in that paper, but again it proved to be quite insensitive to 
configuration, so we will use this value. This introduces an uncertainty in our final value 
for the erasure rate of perhaps ±1 in the exponential.) We can then perform the integral 
in Eq. ([4*4]) numerically. We find that almost all of the erasure occurs in a range of T 
less than 0.5% from the reheat temperature, essentially because v changes very rapidly 
with temperature immediately after the transition. This narrows the available time for the 
Sphaleron erasure, but even for mjj =30 GeV we find that the baryon number is depleted 
by a factor of about exp(11.3). If the conjecture about the effect of Landau damping on 
the negative frequency mode in Q is correct then the suppression is smaller by about 2 in 
the exponent. The results for several Higgs masses are presented in the table; in all cases 
the erasure is very considerable. Since the most optimistic estimates of baryogenesis in the 
minimal standard model can barely account for the current abundance |2(J , this apparently 
rules out electroweak physics as the source of baryogenesis in the minimal standard model. 



5 Conclusion 

It appears that a perturbative treatment of fermions is very well justified near the phase 
transition temperature, and that dimensional reduction should be accurate, though the cor- 
rection from the dimension 6 contribution to the effective potential may not be completely 
negligible. 

How should we understand the results of DPSSG in light of this conclusion? 

First recall what we have done here. We find that at a general \i the fermions will 
induce nonzero corrections to couplings and wave functions, and in particular the correction 
to A is potentially large. Following the idea of the dimensional reduction program we 
combine these contributions with the tree couplings, resulting in renormalization point 
independent couplings At, which are used to compute the Sphaleron energy. The Sphaleron 
energy then already contains that part of the fermion fluctuation determinant which is 
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understood as coupling and wave function corrections; the residual correction, which comes 
from dimension 6 (and higher) operators, is explicitly found to be very small. However, 
had we used the tree couplings at some scale we would then expect fermions to give a 
(/x dependent) nonzero correction due to the difference At — A(/i). Only at one particular, 
convenient renormalization point (Eq. fl42|)) would this contribution vanish. 

Diakonov et. al. use a fixed \i (in |T| they use \i = mw m the proper time scheme, 
and in |2l| they use /i ~ m t exp( r y E /2)) and compute the Sphaleron configuration from the 
tree Lagrangian (corrected however by the thermal contributions to the term). We 
then expect from our work that they should find a correction arising from the difference 
At — A(/x) equal to 

Aitu At - A(//) r°° ,^ 2 h 4 - 1 



9w 9w 



dte—^, (45) 



which will depend on the top mass as gf, exactly as they find. This is not a contradiction 
of the dimensional reduction scheme, which predicts that, because of their use of A(/x) with 
their choice of fi, they should find such a term. We should note that, when this term is 
large (as it is for the choice of \i made in [pl|), one should not trust the functions h, f 
computed from the tree Lagrangian but should include this radiative correction in their 
computation, as we have done here. In other words, the numerical work of Diakonov et. 
al. is probably accurate, but because of the way they have done the problem they will not 
necessarily produce the free energy of the configuration which actually limits the baryon 
erasure rate. 



We should also note that in [II| Diakonov et. al. use a tree, rather than 1 loop, relation 
between m# and A(/i), and a tree, rather than 1 loop, value of A in calculating the phase 
transition temperature. (In |?T) the relation between m H and A(/i) is computed at one 
loop, but the transition temperature is still found using A(/x) rather than At-) This is 
inconsistent with a one loop analysis of the Sphaleron rate and may explain the difference 
in our results for the dissipation of baryon number. 

What is the overall effect of fermions? For small Higgs mass, At > Atr ee) and as 
z/(T reheat ), and hence E Sp h, fall with increasing At, we find more baryon number dissipation 
than we would if we ignored fermions altogether and used A tre e- This is the reason that, 
even for very small Higgs mass, we still find substantial baryon number erasure. We 
should emphasize once more that to find this result it was important to apply one loop 
corrections systematically, in the effective potential, the phase transition temperature, and 
the Sphaleron energy, but that the effect is basically perturbative, and the high temperature 
expansion accounts for it successfully. 

Do the results of the last section preclude electroweak baryogenesis? They make it 
unlikely that baryogenesis can be viable in the minimal standard model. However we 
should note that we have only used the 1 loop effective potential, and while extending our 
results to the two loop potential, which is known, gives essentially the same conclusions, it 
is not clear that the perturbative treatment of the effective potential is reliable; the phase 
transition may be stronger than perturbation theory suggests. Also, we have said nothing 
about extensions to the standard model. For instance, in the two doublet model, the phase 
transition can be much stronger without contradicting experimental bounds on the Higgs 
22], and the Sphaleron bound only narrows the parameter space. 



mass 
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6 Appendix A: Proper Time Calculations 

In this appendix we will show how to compute the finite temperature fluctuation determi- 
nant using a proper time regulation analogous to that of DPSSG. The basic idea is to use 
the relationship 

In = lim -Tr f°° -(e~ Kt - e~ Kot ) , (46) 



DetK Je t 

which holds when K is a quadratic, positive definite operator and K is its free, massless 
approximation. In our case the operator j^d^ — m, or Eq. (|l|), is not quadratic or positive 
definite. Fortunately it has the same spectrum as the operator — 7^<9 M — m, so we can write 
the determinant we actually want in terms of a quadratic, positive definite operator, 

-lnDet^^-m) = — (lnDet^fy - m) + lnDet(-7 M <9 M - m)) 

= -^Tr ln(-7V0,A - Y(d^m) + m 2 ) , (47) 

which acts on a Euclidean space where the time direction is cylcic with antiperiodic bound- 
ary conditions and period T. It therefore already includes all thermal effects, which do not 
have to be added by hand, as in the treatment of DPSSG. However, we will have to be 
careful when we regulate to see that we are making T independent subtractions. 



Following |23] we perform the trace by inserting a complete set of plane-wave states, 



Tie~ Kt = tr / d A xTY, I -^-e- ip x e~ Kt e ip x , (48) 

Po \ / 

where tr is over Dirac indicies. The factor exp(ip-x) can be brought through the operator to 
cancel exp(— ip-x), but in doing this all derivative operators in K are shifted, — > d^ + ip^. 
The zero temperature limit is found by replacing the sum on p with the integral / dpo/(2^). 

We illustrate the renormalization procedure by computing the effective potential in this 
regulation. In this case K = —d 2 + m 2 because the mass is space independent. After the 
shift, the derivatives have nothing on which to act and do not contribute; we can drop 
them. The problem is then to compute 

2 r t i d3x ^ / ^k e ' p2te ~ p2 ° t{e ~ m2t - i} ' (49) 



where we have performed the Dirac trace and removed the trivial integral over dx^ so that 
the result will be the free energy density. 
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Next we expand exp(— m 2 i) in powers of t. The first term is cancelled by the —1. The 
second gives 



- 2 / m 2 d 3 x r dtT V e-fo* / -^-e^ 2 * . (50) 
J Je ^ J {2ny 

This expression is small t divergent, cut off by e. To render the theory cutoff independent 
we should add and subtract a temperature independent expression with the same small t 
behavior and absorb the one we added with a counterterm in the tree level mass. The 
correct expression to subtract is 

dp _ v 2 t f d 3 p 



2/m V,p t /^/^e- , 51 ) 



(2kY 

There is now no obstacle to performing the integral over t, or to setting e equal to zero. 
The result is 



2 W(&( r ?-/t)^ 



Both the sum and the integral over p are straightforward, giving 

1 _ tanhj^ r dp 1 _ 1' 
^ pl+p 2 " 2p ' J 2tt pi + p 2 2p " 1 ' 

combining them and performing the integral over angles, we get 

r 1 r — 1 r m 2 T 2 

_ 2 y* OT »_y__^ = y A _, (54) 

which is the well known expression for the thermal contribution to the Higgs mass squared. 

In their discussion of the renormalization of the theory DPSSG advocate cutting off the 
proper time integral of the counterterm at some finite upper bound, call it ji~ 2 (in their 
case, m^)- Doing so changes the result of the above calculation to 

m 2 T 2 „ 2 m 2 

~12~~^ 2 ~' (55) 

which means there is a discrepancy between the tree level mass squared parameter and 
the renormalized, vacuum mass squared parameter. There is nothing wrong with this in 
principle as long as we remember it is there; we should use the renormalized mass squared 
when performing calculations such as the Sphaleron configuration and remember that we 
have already included part of the fermion contribution by doing so; we will need to subtract 
it off, along with the thermal mass squared. It is much easier and more straightforward, 
however, to follow our procedure and absorb the entire vacuum correction in a counterterm. 
Continuing to expand exp(— m 2 t), the next term is 

+ / m 4 d 3 x rtdtTe-rt [ -^e~ p2 \ (56) 

j Je po j (2ity 

which is logarithmically divergent at small t. Because the lowest Matsubara frequency 
is nonzero, it is cut off exponentially at large t and there are no infrared problems in 
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its evaluation. Again we have an available Lagrangian parameter in which to absorb the 
divergence, and we should add and subtract a T independent expression with the same 
ultraviolet behavior, and absorb the one with the same sign into the Higgs self-coupling 
parameter. The corresponding vacuum integral is both ultraviolet and infrared divergent, 
and to prevent infrared divergences we must introduce a renormalization scale into the 
problem. Following DPSSG, we subtract 



rd Po »h f d3 P 



giving 



m 4 d 3 x / tdt / ^e" p o* / j^e- p 1 , (57) 
2n J (2n) 6 



jr d 3 p give.' 

and gives + ln(/i 2 e)/(167r 2 ). The remaining integral, 



(27T) 

The integral over d 3 p gives t~ 3 l 2 /(8tc 3 ^ 2 ). The counterterm can then be evaluated directly 



T 

57T2 po 



t-2dte- p °\ (59) 



is more delicate. For small terms in the sum the integral over t is approximately ^/tt/po, 
so the early terms in the sum are 

T 1 

(60) 



47r j=1X i7rT 



The difference between this sum and a corresponding integral is concentrated in the first 
few terms. The sum is approximately 



1 

47T 



dpo 1 \ 7b + In 2 



(61) 



/ttT 27T Po / 27T 

so smoothing over the sum in Eq. (|59"D introduces a correction of (7^ + ln2)/(87r 2 ), giving 

8vr 2 + 4vr f It 2vrp 7 £ ' l ° j 

The integral over po can be performed by parts, giving 

1E + + -4^ (W + i Ine - 4(1/2)) , (63) 



8tt 2 4^1 2^ V 2 2 

with ip the digamma function, i[)(l/2) = — 7^ — 21n2. The result is then 

m4d3x7E -21n, + ln(^) (M) 

The fraction in the integral gives the proper time renormalization value of DA. 

The terms higher order in m 2 are actually easier; the term at order m 2n is given by the 
integral 



* 7 n! (2tt) 3 v ; 
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Performing the integral over t, we get 



n 

The result per unit volume is 

(-l) n m 2n 



2("l) n f ,3 2n^ f°°P 2 dp 1 , RR , 

a xm > / — r-s or— • 66 

^ 'o 2tt 2 (p 2 + pl) n v ; 




(67) 



from which Eq. (J|) follows immediately by performing the sum and the integral. Note that 
the calculation was completely large t finite because the exp(— p^t) term always decays 
exponentially fast, since pg is always at least (7rT) 2 . This is to be contrasted with the zero 
temperature [^3 and zero Matsubara frequency bosonic [15] cases and explains why the 



expansion in operator dimension is possible at finite temperatures for fermions, while it is 
known to have trouble at zero temperature and for the zero Matsubara frequency bosonic 
modes. 

The same basic techniques can be used for the case with gauge fields and spatial vari- 
ation. A complication which arises when computing other dimension 4 operators is that, 
in contributions involving A , powers of po appear in addition to exp(— p^t). The result is 
that the argument of the digamma function shifts by 1 per power of Pq. This is why the 
dimension 4 terms containing Aq are not simply multiples of D4. 

Next we will compute the fermionic contribution to the vacuum effective potential in 
this regulation. From the above, the contribution from each color is 



dt f d* P e - v H {e -mH _ 1 + mH) ) _ ( f 2 dt i 



t J (2tt) 4 v 7 \Je t J (2tt) 4 2 



(68) 



The —1 is the vacuum energy subtraction, the m 2 t is the mass squared counterterm, and 
the last expression is the self-coupling counterterm. The integrals over p may be performed 
immediately, giving 




, \ 1 m 4 ft* 2 dt 



After integrating the first expression by parts three times, we get 



m 4 / 3 , a 2 \ 

1 + In ^- - lE . (70) 



167T 2 \2 m? 

Recalling that m 2 = g 2 &Q, summing on colors, and adding this expression to the tree level 
effective potential 

- m2$t$ + A($ f $) 2 (71) 

we find that the curvature of the effective potential (the second derivative with respect to 
v) at the minimum [y = v a ) is 

V" = 2u 2 (x + M(\n^- lE )) , (72) 
V 167r^ mf I 
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and the effective potential parameter TOq is 



2 v N c9 y 
< = T + ^- (73) 

The analogous MS result is the same but with the — je removed from Eq. fl72p. 

V" is not the physical Higgs mass squared. The physical Higgs mass is the ratio of 
potential to kinetic energy of a long wavelength fluctuation about the minimum, times 
the wave number of the fluctuation; there is a wave function correction which requires 
computing the self-energy at the pole mass. We have been unable to do the calculation in 
proper time regulation (see however pl|), but we have performed it in MS using an integral 
from p4| . The result is that 



The proper time value is presumably the same but with a — ^ E inserted after the log. For 
4m 2 >> m 2 H the expression following the log is about 2m 2 H /3(4:m 2 — m 2 H ) ~ and we have 
permitted ourselves to neglect it when relating the quartic coupling to the physical Higgs 
mass. 



7 Appendix B: Other high dimension operators 

We list here those dimension 4 and 6 operators induced by fermions which contain both weak 
577(2) and strong S77(3) fields. These are probably of no phenomenological consequence 
and almost certainly do not influence the strength of the phase transition. We include them 
only for completeness. 

We write the time component of the gluon field as A Qg and the field strength tensor as 
Gij. Denoting the number of families as Ne = 3, the dimension four, mixed terms are 



9t9 s A 2 flst^ 9 w 9 s Nf a 2 a2 



(75) 



The dimension six mixed terms are 

,2^2 „2 

,2„2 AT_na „2„2 



There are also pure glue terms which are identical to the pure gauge terms listed in the 
text, but with the substitution —>■ 2Nf where Nf = 6 is the number of quark flavors. 
The 2 is because both right and left handed quarks couple to the gluons, wereas the weak 
coupling is chiral. 
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m H (GeV) 


30 


50 


60 


Trchcat(GeV) 


71.2127 


86.6495 


95.6517 


Esph/T 


29.18 


24.79 


22.38 


E Sph /T from 


($t$)3 


1.6 


0.60 


0.31 


$T$(A$) T A$ 


-6.3 x 10~ 3 


-3.7 x 10~ 3 


-2.7 x 10~ 3 




-3.3 x 10~ 3 


-1.8 x 10~ 3 


-1.3 x 10~ 3 


(A$) T $$ T A$ 


-2.1 x 10~ 3 


-1.2 x 10~ 3 


-8.6 x 10" 4 


F. a F. a $t(T) 


-9.9 x 10~ 4 


-6.2 x 10~ 4 


-4.7 x 10" 4 


(D i F ij ) a (D k F kj ) a 


-4.6 x 10~ 4 


-2.9 x 10" 4 


-2.2 x 10~ 4 


f pa pb pc 
J abc r ij r jk r ki 


-8.4 x 1(T 5 


-5.3 x 10~ 5 


-4.0 x 10~ 5 


(D 2 ®)W 2 <S> 


-4.1 x 10"° 


-4.2 x 10- y 


-4.4 x 10 _e 


deplel 


Aon factor 


exp(11.3) 


exp(15.0) 


exp(17.1) 


m H (GeV) 


70 


80 


90 


^reheat (GeV) 


105.1149 


114.8275 


124.6283 


Esph/T 


20.49 


19.26 


18.70 


Esph/T from 


($t$)3 


0.17 


0.100 


0.065 




-2.0 x 10~ 3 


-1.6 x 10" 3 


-1.5 x 10~ 3 




-8.8 x 10~ 4 


-6.7 x 10~ 4 


-5.6 x 10~ 4 


(A$) T $$ T A$ 


-6.3 x 10~ 4 


-4.9 x 10~ 4 


-4.3 x 10~ 4 


F. a F. a $t(T) 


-3.6 x 10" 4 


-3.1 x 10" 4 


-2.9 x 10" 4 


(D t F tj ) a (D k F kj ) a 


-1.7 x 10~ 4 


-1.5 x 10~ 4 


-1.4 x 10" 4 


f jpa rpb rpc 
J abc r ij r j k r k i 


-3.1 x lO" 5 


-2.6 x 10~ 5 


-2.5 x 10~ 5 


{D 2 <&)W 2 <$> 


-4.8 x lO" 


-5.5 x 10~ d 


-7.0 x 10-° 


deplel 


ion factor 


exp(19.0) 


exp(20.3) 


exp(21.5) 



Table 1: Reheat temperature, Sphaleron energy at the reheat temperature, and contribution 
of dimension 6 operators to Sphaleron energy for a number of physical Higgs masses 



Figure 1: Feynman diagrams corresponding to two insertions of Hj. 
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